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ABSTRACT 

We use the GEANT4 Monte Carlo framework to calculate the 7-ray albedo of the Moon due to interactions of 
cosmic ray (CR) nuclei with moon rock. Our calculation of the albedo spectrum agrees with the EGRET data. 
We show that the spectrum of 7-rays from the Moon is very steep with an effective cutoff around 3^1 GeV 
(600 MeV for the inner part of the Moon disk) and exhibits a narrow pion-decay line at 67.5 MeV, perhaps 
unique in astrophysics. Apart from other astrophysical sources, the albedo spectrum of the Moon is well 
understood, including its absolute normalisation; this makes it a useful "standard candle" for 7-ray telescopes. 
The steep albedo spectrum also provides a unique opportunity for energy calibration of 7-ray telescopes, such 
as the forthcoming Gamma Ray Large Area Space Telescope (GLAST). Since the albedo flux depends on the 
incident CR spectrum which changes over the solar cycle, it is possible to monitor the CR spectrum using the 
albedo 7-ray flux. Simultaneous measurements of CR proton and helium spectra by the Payload for Antimatter- 
Matter Exploration and Light-nuclei Astrophysics (PAMELA), and observations of the albedo 7-rays by the 
GLAST Large Area Telescope (LAT), can be used to test the model predictions and will enable the LAT to 
monitor the CR spectrum near the Earth beyond the lifetime of the PAMELA. 

Subject headings: elementary particles — line: formation — radiation mechanisms: non-thermal — Moon — 
cosmic rays — gamma-rays: theory 



1. INTRODUCTION 

Interactions of Galactic CR nuclei with the atmospheres of 
the Earth and the Sun produce albedo 7-rays d ue to the de- 
cay of secondary neutral p ions and kaons (e.g., ISeckel et al.l 
1 199 It lOrlando et all 120071) . Similarly, the Moon emits 7- 
rays due to CR intera ctions with its surface (Morris 1984; 
ITho mpson et al. 1997). However, contrary to the CR inter- 
action with the gaseous atmospheres of the Earth and the 
Sun, the Moon surface is solid, consisting of rock, making 
its albedo spectrum unique. 

Due to the kinematics of the collision, the secondary par- 
ticle cascade from CR particles hitting the Moon surface at 
small zenith angles develops deep into the rock making it dif- 
ficult for 7-rays to get out. The spectrum of the albedo 7-rays 
is thus necessarily soft as it is produced by a small fraction of 
low-energy splash particles in the surface layer of the moon 
rock. The high energy 7-rays can be produced by CR par- 
ticles hitting the Moon surface with a more tangential trajec- 
tory. However, since it is a solid target, only the very thin limb 
contributes to the high energy emission. 

The 7-ray albedo of the Moon has been calculated by 
Morris (1984) using a Monte Carlo code for cascade de- 
velopment in the Earth's atmosphere that was modified for 
the Moon conditions. However, the CR spectra used as in- 
put in the iMorrisI calculat ion differ considerab ly from recent 
measurements by A MS (lAlcaraz et all l2000allbh and BESS 
(ISanuki et al.l 1200 0) at both low and hi gh ener gies. In par- 
ticular, the CR proton spectrum used bv lMorrisI produces the 
correct intensity only around 10 GeV. Additionally, due to 
the lack of accelerator data and models a number of approx- 
imations and ad-hoc assumptions were required to calculate 

1 Also Kavli Institute for Particle Astrophysics and Cosmology, Stanford 
University, Stanford, CA 94309 



the hadronic cascade development in the solid target of the 
Moon's surface. 

The Moon has been detected by the EGRET as a point 
source with integral flux F(>100 M eV) = (4.7±0.7) x 10~ 7 
cm~ 2 s " 1 (IThompson et"alj|199 7). ~24% below the predic- 
tions bv lMorrisI (1984) although the spectral shape agrees with 
the data. The observed spectrum is steep and yields only 
the upper limit ~5.7 x 10~ 12 cm~ 2 s _1 above 1 GeV. At 
lower energies, 7-ray spectroscopy data acquired by the Lunar 
Prospector have been used to map the elemental composition 
of the Moon surface (lLawrence et alJlT998h iPrettyman et all 
120061) . 

In this paper, we calculate the 7-ray albedo from the Moon 
using the GEANT4 software framework and discuss the con- 
sequ ences of its measurement by the upcoming GLAST mis- 
sion ([Michelson 2007; Ritz 20 07|). Our pr eliminary results are 
presented in Moskalenko & Porterj (120071) . 

2. MONTE CARLO SIMULATIONS 

Sin ce the first calculation of the Moon albedo (IMorrisI 
119841) . computers and Monte Carlo codes have advanced 
considerably, while more data on particle interactions and 
CR flux es have become avail able. The modern GEANT4 
toolkit 2 ( Agostinelli et al. 2003) is widely used across par- 
ticle physics, medical science, and space physics to simu- 
late the passage of particles through matter. It provides ex- 
tensive functionality including physics models, tracking, ge- 
ometry, and visualisation. A full range of electromagnetic 
and hadronic physics processes over a wide energy range are 
available including long-lived particles. A large set of mate- 
rials and elements are also included. Together with the ex- 
tensive documentation and examples, it is straightforward to 

2 Available from http://cern.ch/geant4 
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FIG. 1. — Beam/target/detector setup for simulating CR interactions in 
moon rock. The primary beam enters the moon rock target with incident polar 
angle 8 P . Secondary 7-rays are emitted with polar angle 8. The detection 
volume surrounds the target. 



describe and simulate complicated detectors and interactions. 

In the present work, we use version 8.2.0 of the GEANT4 
toolkit. Figure [T]illustrates our beam/target/detector setup for 
simulating CR interactions in the Moon. The primary CR 
beam (protons, helium nuclei) is injected at different incident 
angles into a moon rock target. We take the composition of the 
moon rock to be 45% Si0 2 , 22% FeO, 1 1 % CaO, 10% A1 2 3 , 
9% MgO, and 3% Ti02 by weight, consistent with mare basalt 
meteorites and Apollo 12 and 15 basalts (|Lawrence et al.l 
1998l lAnand et aU2003l:|Prettvman et ail2~0 06). A thin hemi- 
spherical detector volume surrounding the target is used to 
record the secondary 7-ray angular and energy distributions 
in the simulation. 

In GEANT4 the specification of the physics list is of consid- 
erable importance since it is where all the physics to be used 
in the simulation is defined. In the physics list, the user de- 
fines all necessary particles in the simulation, assigns physics 
processes to each, and sets appropriate range cuts for sec- 
ondary production. In GEANT4 there is a distinction between 
processes and models: processes take care of particle decays 
and interactions while for any process there may be one or 
more models and cross section parameterizations. For most 
processes (electromagnetic, elastic, decay) there is only one 
model. However, for hadronic processes and, in particular, 
inelastic processes, more than one model can be used over 
different energy ranges. The preference for one model over 
another is not specified in the documentation but depends on 
the particular application. Fortunately, to guide users of the 
toolkit, extensive examples are provided for use cases in dif- 
ferent settings. 

To construct our physics list, we examined the two exam- 
ples in the GEANT4 distribution that most closely resembled 
our simulation setup: the "cosmicray .charging" and "radio- 
protection" examples in the 8.2.0 release. We included all rel- 
evant particles (leptons, photons, pions, kaons, heavy mesons, 
protons, antiprotons, neutrons, antineutrons, and resonances) 
and included the standard electromagnetic processes (ionisa- 
tion, bremsstrahlung, multiple scattering, etc.) appropriate for 
each particle. For hadronic processes, we included "at rest" 
processes, such as negative pion and kaon absorption, neutron 
and muon capture, antiproton and antineutron annihilation, 
and "in flight" processes such as elastic and inelastic scatter- 
ing (photo- and electro-nuclear and other secondary produc- 
tion), in-flight capture, and fission. For all long-lived unstable 
particles the decay process was also included. 



For inelastic processes GEANT4 provides different types 
of hadronic shower models, each of which has an intrinsic 
range of applicability. Three types of hadronic shower models 
are available: data-driven models for high precision treatment 
of the transport of low energy neutrons, parameterised mod- 
els based on a re-implementation of the GHEISHA package 
of GEANT3.21, and theory -based models that cover a wide 
energy range. Pre-packaged combinations of these models 
are provided with the GEANT4 distribution. However, for all 
of these inelastic scattering by ions is only treated up to 100 
MeV. Cosmic-ray helium is a significant component of the CR 
flux interacting in the Moon's surface, and we consider beam 
energies considerably higher than this in our simulations. Fur- 
thermore, many of the pre-packaged lists are intended for use 
in applications such as high energy collider simulations where 
the treatment of low energy processes is less important. 

Guided by the aforementioned examples and the documen- 
tation, we constructed the hadronic physics list using the fol- 
lowing models. For energies up to 70 MeV, a pre-equilibrium 
decay and evaporation (PreCo) model for interactions of pro- 
tons and neutrons was used. Above this, a Bertini cascade 
(BERT) model between 60 MeV and 4 GeV, a parameterised 
model (LEP) from 3 GeV to 15 GeV, and a quark gluon string 
(QGS) model from 12 GeV up to 100 GeV were employed. 
For pions and kaons, the BERT model was used up to 4 GeV 
with the LEP and QGS models used over the same energy 
ranges as for protons and neutrons. Deuterons, tritons, and al- 
pha particles use the low energy inelastic code up to 100 MeV, 
and a binary ion cascade model from 40 MeV to 40 GeV. Each 
model is used in its applicable energy range (D. Wright, pri- 
vate communication). In the energy ranges where there are 
overlaps, the GEANT4 code interpolates between the mod- 
els. The overlap ranges were chosen to ensure as smooth 
a transition between the models as possible. These models 
compare well with data and other particle tra nsport codes in 
standard accelerator and space environments (Wellischl l2003l: 
lErsmarketai1l2004l) 3 . 

The 7-ray yield dYy(E p , cos 9 p ) / dE^d cos 9 is calculated 
using the GEANT4 beam/target setup with a Monte Carlo 
method; here E p is the kinetic energy per nucleon of the inci- 
dent particle, 9 P is the incident polar angle, E 1 is the energy 
of 7-rays, and 9 is the polar photon emission angle. 

Figure|2](left) shows the total secondary 7-ray yield and the 
contribution by different processes integrated over all emis- 
sion angles outward from the Moon's surface for protons with 
E p = 10000 MeV and incident angle cos6> p = 0.1. The 7- 
ray emission is produced in a number of processes: pion and 
kaon decay, secondary electron and positron bremsstrahlung, 
and so forth. A considerable flux of 7-rays is produced in nu- 
clea r reactions such as neutron capture and none lastic scatter- 
ing (lLawrence et al.ll998l:lPrettyman et al.l2006l) : the features 
below ^10 MeV are due to nuclear de-excitation lines, where 
the most prominent contribution comes from oxygen nuclei. 
The development of the cascade in the rock causes modifica- 
tion of the yield spectra from the usual spectral distributions. 
For example, the 7r°-decay spectrum would normally be sym- 
metric about 67.5 MeV, instead it has a significant low-energy 
extension due to Compton scattering of the 7-rays in the rock. 

Figure|2](right) shows the secondary 7-ray yields integrated 
over all emission angles outward from the Moon surface for 
protons with E p = 500 MeV and 10000 MeV at incident an- 

3 see also physics lists at http://geant4.slac.stanford.edu/ and 
http://www.particle.kth.se/desire/ 
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FIG. 2. — 7-ray yield per proton interaction integrated over all emission angles from the Moon surface. Left panel: yield calculated for E p = 10000 MeV 
and cos 8 p = 0.1 with components shown. Line styles: blue-dashed, decay; red-long-dash-dot, bremsstrahlung; magenta-short-dash-dot, non-elastic scattering 
and de-excitation; cyan-dot, low energy neutron capture; black-solid, total. For each component, the initial process responsible for the production of the 7-ray 
secondary is given by the line style but the yield distribution includes also processes such as Compton scattering. Right panel: yield calculated for two different 
energies and incident angles. Line styles: blue-solid, cos 9 P = 0.1; red-dashed, cos 9 P = 1. Line-sets: lower, E p = 500 MeV; upper, E p = 10000 MeV. 



gles cos P = 0.1 and 1, respectively. The shapes of the yield 
curves for different incident angles are very similar to each 
other for the case of low energy protons where the secondary 
particles (pions, kaons, neutrons, etc.) are produced nearly at 
rest. In the high energy case, the secondary distribution for 
protons incident near zenith has a cutoff above ^500 MeV. 
Further away from zenith, higher yields of secondary 7-rays 
are produced while the spectrum of 7-rays becomes progres- 
sively harder. This is a result of the cascade developing mostly 
in the forward direction: for near zenith primaries, most high 
energy secondary 7-rays will be absorbed in the target, while 
a small fraction of splash albedo particles, mostly low energy 
ones, produce the soft spectrum 7-rays; further from zenith, 
the high energy secondary 7-rays will shower out of the Moon 
surface. 

The clear narrow line at 67.5 MeV in Figure [2] (left and 
right panels) is due to neutral pions decaying near rest. 
To ensure that it is not an artefact of the simulation, we 
made an extensive analysis of the hadronic interaction mod- 
els employed in GEANT4 and various geometries of the 
beam/target/detector configuration. Independently of the 
hadronic interaction model used, the line is always present, 
but its intensity may vary slightly. However, dependent on 
the beam/target/detector configuration the line may appear or 
disappear. The most feasible explanation of this effect is the 



following. The line is always present since there is always 
some small fraction of pions decaying at rest, but in the most 
common configurations, such as thin target, gaseous target, 
and/or detector position in the forward hemisphere, the line is 
hidden under the strong continuum background. In our case of 
a solid thick target and the detector location in the backward 
hemisphere, a considerable part of the continuum background 
is cut off and the line becomes visible. 

There are several reasons why the continuum background 
is cut off. The interaction length of a GeV proton in oxy- 
gen (the most abundant element in the moon rock) is ^80 g 
cm~ 2 , which corresponds to ^25 cm for the moon rock (den- 
sity ~3.3 g cm" 3 ). As the cascade develops, energetic parti- 
cles penetrate deep into the rock, too deep for 7-rays to come 
out. The 7-ray albedo is produced by relatively low-energy 
particles at the depth of about one radiation length, ~34 g 
cm~ 2 , which corresponds to ^10 cm for the moon rock. The 
pions decaying in the upper layer of the moon rock are nec- 
essarily low-energy (splash) particles while the high-energy 
part of the cascade producing the continuum spectrum is cut 
off because of the kinematics of the interaction and the for- 
ward cascade development. Since many photons are coming 
from a large depth, this explains the large contribution of the 
Compton scattered 7-rays below ^60 MeV (Figure |2). For 
lower energy protons, this also explains why the 67.5 MeV 
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line is not prominent, since the development of the cascade 
then occurs relatively close to the surface producing a contin- 
uum that essentially washes out any line feature. 

3. CALCULATIONS 

The CR spectrum above the geomagnetic cutoff near Earth 
(at 1 AU) can be directly measured by balloon-borne instru- 
ments or spacecraft. However, it has been done during short 
flights at different phases of solar activity. To calculate the 
Moon albedo at an arbitrary modulation level, we use the 
local interstellar (LIS) spectra of CR protons and helium as 
fitted to the numerical results of the GALPROP propagation 
m odel (reacceleration and plain diffusion models, Table 1 
in iPtuskin et al.1 (120061) '): the CR particle flux at an arbitrary 
phase of solar activity at 1 A U can then be estimated using 
the force-field approximation dGleeson & Axfordll 19681) : 



dJ p (E k ) _ dJ™(E k +$Z/A) 



E 2 - M 2 



dE k 



dE k 



{E + <S>Z/A) 2 - M 2 ' 



where dJ^° / dEk is the LIS spectrum of CR species, E k is the 
kinetic energy per nucleon, E is the total energy per nucleon, 
$ is the modulation potential, Z and A are the nucleus charge 
and atomic number correspondingly, and M is the nucleon 
mass. 

To fit the LIS CR spectra we choose a function of the form: 



dJ c 



dE k 



Ja^a^Ek+h)- 



(2) 



where the flux units are m~ 2 s _1 sr _1 (GeV/nucleon) -1 and 
the parameter values are given in Table [TJ The latter are 
not unique and other sets could produce similar quality fits, 
but this does not affect the final results. Figure [3] shows 
the LIS CR pro ton spectrum and modulated one compared 
with th e BESS dSanuki et alj |2000) and AMS dAlcaraz et alj 
2000a) data taken during the period of moderate solar activ- 
ity. The parameters given in Table Q]yield a LIS proton spec- 
trum that is somewhat higher than the GALPR OP LIS be- 
low ~ 2 GeV (plain diffusion model 44,999726. IPtuskin et ail 
120061) . However, the fit to the AMS and BESS data after 
modulation ($ = 550 MV) is good. The spectrum labelled 
$ = 1500 MV corresponds to the period of high solar ac- 
tivity. The proton spectrum used by Morris ( 1984) is shown 
by the thick solid line. The disagreement with the modern 
CR proton data is considerable. For the helium spectrum, the 
parameters given in Table Q] provide a good fit to the GAL- 
PROP LIS, and agree with the AMS dAlcaraz et al.ll2000bl) 
data after modulation ($ = 550 MV); at high ener gies the fit 
also passes through JA CEE (|Asakimori et al. 1998) and Sokol 
dlvanenko. et al.lll993l) data. 
The 7-ray albedo flux at the Earth is calculated as: 
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where ip is the angular distance from the geometrical center of 
the Moon disk as seen from the Earth, 1 / cos is the Jacobian, 
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FIG. 3. — The CR proton spectrum. The dashed lines are the LIS spectrum 
(upper) and modulated one (lower, <E> = 550 MV) as calcula ted by GAL- 
PROP in the plain diffusion model 44.999726 JPtuskin et all2006ft . The thin 
solid lines are the fit to the LIS spectrum (upper) and modulated spectra (mid- 
dle, 5> = 550 MV; low er <S> = 1500 MV). The thick solid li ne is the CR 
proto n spectrum used by Morris 1 1984 ). Data: circles — BESS I San uki et all 

" "Ha - 



l2007j) , triangles — AMS I Alcar azetaU2000al) 



R = 1738.2 km is the Moon radius, and r = 384401 km is 
the distance to the Moon. 

Figure [4] shows the calculated total 7-ray albedo spectrum 
for CR protons and helium compared to the EGRET data for 
periods of lower (upper solid: $ = 500 MV) and higher 
(lower solid: $ = 1500 MV) solar activity. Taking into ac- 
count that the exact CR spectra during the EGRET observa- 
tions are unknown, the agreement with the data is remarkable. 
The broken lines show the spectra from the limb (outer 5') and 
the central part of the disk (20' across) for the case of higher 
solar activity. As expected, the spectra from the limb and the 
central part are similar at lower energies (<10 MeV); at high 
energies the central part exhibits a softer spectrum so that vir- 
tually all photons above ^600 MeV are emitted by the limb. 

Interestingly, the pion-decay line is still significant even 
when the 7-ray yields are integrated over the spectrum of CR, 
although, its intensity is small relative to the total albedo flux. 
The line perhaps always exists due to the splash albedo pi- 
ons, but in case of a thin and/or gaseous target (usual in as- 
trophysics) it is indistinguishable from the background 7-rays 
(see Section|2]for more discussion). 

4. DISCUSSION 

The GLAST LAT is scheduled for launch by NASA in early 
2008. It will have superior angular resolution and effective 
area, and its field of vie w (FoV) will far exce ed that of its pre- 
decessor, the EGRET dMcEnerv et alj|2004l) . The LAT will 
scan the sky continuously providing complete sky coverage 
every two orbits (approximately 3 hr). The on-axis effective 
area of the LAT increases from ^3000 cm 2 at 100 MeV to 
-8500-9000 cm 2 at 1 GeV and higher 4 . In this case, the 
point spread function (PSF) of the instrument has a 68% con- 
tainment radius —4.3° at 100 MeV reducing dramatically at 
higher energies: -0.8° at 1 GeV, -0.5° at 2 GeV, and -0.2° 
at 10 GeV. Using only events from the front section of the 
LAT, the PSF improves to a 68% containment radius of ~3° 
at 100 MeV, -0.5° at 1 GeV, ~0.25° at 2 GeV, and ~0.1° at 



see http://glast.stanford.edu 
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Table 1 . Fits to local interstellar CR spectra 
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NOTE. — Intensity units: m 2 s 1 sr 1 (GeV/nucleon) 1 . 
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FIG. 4. — Calculated 7-ray albedo spectrum of the Moon. Line-styles: 
black-solid, total; blue-dotted, limb - outer 5'; red-dashed, center — inner 
20'. Upper solid line: $ = 500 MV; lower solid line: $ = 1500 MV. 
Limb and center co mponents are only s hown for <E> = 1500 MV. Data points 
from the EGRET I Thompson et al. 1997) with upper and lower symbols cor- 
responding to periods of lower and higher solar activity, respectively. The 
differential 1 year sensitivity of the LAT is shown as the shaded region. 



10 GeV, but the on-axis effective area is essentially halved. 

About 20% of the time the Moon will be in the FoV at 
different viewing angles. The photon flux expected from the 
Moon for $ = 500 MV above 100 MeV, 1 GeV, and 4 GeV is 
-5 x 10~ 7 , -2 x 10~ 8 , and <1.6 x lO" 10 photons cm~ 2 s _1 , 
respectively. For $ = 1500 MV, the expected flux is reduced 
by a factor of 2 at 100 MeV only. With these fluxes and the 
above values for the effective area, and allowing for an addi- 
tional factor 2 reduction to take into account time off-axis, in- 
strumental deadtime, and South Atlantic Anomaly traversals, 
we estimate the LAT will collect ~5 x 10 3 , ~6 x 10 2 , and <5 
photons, respectively, in the above energy ranges in one year. 
These numbers are reduced by a factor of two if only events 
from the front section of the LAT are used. Interestingly, the 



albedo flux at low energies is high enough that the impact of 
the broader PSF is significantly reduced. 

Measuring 7-rays from the Moon presents several interest- 
ing possibilities for the LAT. Apart from other astrophysical 
sources, the albedo spectrum of the Moon is well understood, 
including its absolute normalisation, while the Moon itself is 
a "moving target" passing through high Galactic latitudes and 
the Galactic centre region. This makes it a useful "standard 
candle" for the GLAST LAT at energies where the PSF is 
comparable to the Moon size 0.5°, i.e., at ~1 GeV and higher. 
At these energies the albedo flux is essentially independent of 
the solar modulation. At lower energies the 7-ray flux de- 
pends on the level of the solar modulation and thus can be 
used to infer the incident CR spectrum. A simultaneous pres- 
ence of the PAMELA on-orbit capa ble of measuring p rotons 
and light nuclei with high precision dPicozza et al. 2006|) pro- 
vides a necessary input for accurate prediction of the albedo 
flux and a possible independent calibration of the GLAST 
LAT. An additional bonus of such a calibration is the possi- 
bility to use the GLAST observations of the Moon to monitor 
the CR spectra near the Earth beyond the projected lifetime of 
the PAMELA (currently 3 years). 

The line feature at 67.5 MeV from 7r°-decay produced by 
CR particles in the solid rock target is interesting. There is 
no other astrophysical object predicted to produce such a nar- 
row line and there is no other line expected except, perhaps, 
from dark matter annihilation. The lower energy limit of the 
LAT instrument is below 20 MeV while the energy resolution 
is ~15% at 100 MeV, and improves at higher energies. With 
a suitable event selection it may be possible to observe the 
line; if so, it will provide a possibility of in-orbit energy cal- 
ibration. Another possibility for energy calibration at higher 
energy is provided by the steep albedo spectrum above 100 
MeV: a small error in the energy determination will result in 
a large error in the intensity. 

Because of the steep albedo spectrum, very few photons 
will be detected above ~3— 4 GeV. The central part of the 
Moon has an even steeper spectrum with an effective cutoff 
at ~600 MeV. Above a few GeV, the Moon presents almost a 
black spot on the 7-ray sky providing an opportunity to screen 
out a piece of the sky. Additionally, the expected irreducible 
background rate 5 at these energies within the solid angle sub- 
tended by the Moon is ~0.1 photons per year (S. Digel, pri- 
vate communication). The absence of any significant excess 
over the predicted rate will be a useful cross check on the ex- 
pected irreducible background. 

5. CONCLUSIONS 

The 7-ray albedo of the Moon makes it a unique calibra- 
tion target for 7-ray telescopes. Apart from other astrophysi- 
cal sources the albedo spectrum is well understood, including 

5 7-rays produced in the FoV by CR interactions in the inert material sur- 
rounding the LAT. 
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its absolute normalisation. The albedo of the Moon is dim, 
especially its central part at high energies. Its lower energy 
part exhibits a narrow pion-decay line at 67.5 MeV, perhaps 
unique in astrophysics and never before observed, while its 
continuum intensity depends on the phase of the solar cycle 
and allows one to monitor the ambient spectrum of CR par- 
ticles. The GLAST LAT instrument is well suited for such 
observations. 
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